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Abstract 

We analyze the problem of calculating the solutions and the spectrum of a 
string with arbitrary density and fixed ends. We build a perturbative scheme 
which uses a basis of WKB-type functions and obtain explicit expressions for 
the eigenvalues and eigenfunctions of the string. Using this approach we show 
that it is possible to derive the asymptotic (high energy) behavior of the string, 
obtaining explicit expressions for the first three coefficients (the first two can 
also be obtained with the WKB method). Finally using an iterative approach we 
also obtain analytical expressions for the low energy behavior of the eigenvalues 
and eigenfunctions of a string with rapidly oscillating density, recovering (in a 
simpler way) results in the literature. 

Keywords: Hclmholtz equation; inhomogeneous string; perturbation theory; 
collocation method 



1. Introduction 

In this paper we consider the problem of calculating the normal modes of a 
string with arbitrary density and finite length. In a recent paper, ref. we have 
considered the same problem, and we have proved the following results: first, 
that it is possible to apply iterative techniques and obtain increasingly refined 
approximations to selected normal modes of the string, starting from a particular 
ansatz; second, that it is possible to use a perturbative approach, corresponding 
to an expansion in the small inhomogcncities of the string, and thus prove 
that to second order in this expansion the perturbative results reproduce the 
asymptotic high energy behavior obtained with the WKB method. This result, 
which is not a-priori obvious, depends crucially on the dimensionality of the 
problem, since in one dimension the perturbative contributions depending on 
off-diagonal terms are not neglegible in the limit of asymptotically large energies 
(as it seems to be the case in more than one dimension). 

The purpose of this work is twofold: on one hand, to devise an alterna- 
tive perturbative approach, which is applicable to strings of arbitrary density 
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(including highly inhomoegeneous cases); on the other hand to show that the 
iterative theorems of ref. [J may be used to obtain extremely precise (numerical 
and analytical) results when applied to non-trivial problems. 

The paper is organized as follows: in Section[3]we review the density pertur- 
bation theory (DPT) of ref. [1| ; in Section[3]we devise an alternative perturbative 
method which uses a basis built out of WKB solutions and derive explicit an- 
alytical expressions for the first three coefficients of the asymptotic expansion 
for the energy: while the first two coefficients were already known (they can 
be obtained with the standard WKB method), the third coefficient is new (to 
the best of our knowledge) and it takes the form of an infinite series of diver- 
gent series (each of which can be Borel resummed); in Section @] we describe a 
modified version of this alternative perturbative method, which uses WKB-type 
solutions, corresponding to non-physical densities and which may provide a bet- 
ter description of selected modes of the string; in Section [S] we briefly review 
the iterative theorems of ref. [l| ; in Section [5] we discuss several applications of 
the methods of previous sections: in particular we test with high precision the 
analytical formula for the third asymptotic coefficient, and we apply the itera- 
tive methods to obtain the asymptotic behavior of the low energy spectrum of 
a string with highly oscillating density, thus reproducing results obtained in the 
literature |2|; finally in Section [5] we draw our conclusions. 



2. Density perturbation theory 

The starting point of our discussion is the Helmholtz equation 



dx 2 



E n p{x)^ n {x) 



(1) 



which describes the vibrations of a string of density p(x). We assume that the 
string is centered in the origin and that it has a length 2L. The solutions \P n (x) 
are orthogonal [3| 



y„(x)ty m (x)p{x)dx = <5„ 



(2) 



and obey specific boundary conditions at the border. In this paper we consider 
only Dirichlct boundary conditions, i.e. \P n (±L) = 0. 

Equation (flj can be cast into a more convenient form by introducing the 
functions $„(x) = p{x)^ n {x) and writing it as 



The operator O = 



d?_ 

dx 2 

l 

fp{x) 



_ d 2 

dx' 2 



<f> n (x) = E n $ n (x) 



(3) 



/p(x) 



has a manifestly hermitian form 



and its eigenvalues and eigenfunctions are directly related to those of eq. ([T]). 
The results of refs.0, 0, 0] are obtained using this symmetrized form for the 
differential operator of the inhomogeneous Helmholtz equation. 
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We now assume that the density p(x) of the string may be expressed as 



p{x) = p {l + rj Sp{x)) 



(4) 



where |<5/c(a;)| <C 1 and p(x) > for x € (-L, L). 77 is a non-physical parameter 
used to allow power-counting (eventually this parameter is set to one). Un- 
der these assumptions dp represents a perturbation and therefore one may use 
perturbation theory. 

We expand the operator O in a power series 



where 0,13] 
PoOq 
P0O1 
PoO 2 
PoO 3 



dx 2 
1 



dp 



25p 



3 

16 



6 P 2 



5 
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n=0 



O n rf 



dx 2 

dx 2 r 

dx 2 
d 2 
dx 2 



6p 



dx 2 
+ 36p 2 

Sp + Sp 



5p 
dx 2 

jL 

dx 2 



dx 2 



dx 2 



SP 2 



SP 1 



(5) 

(6) 
(7) 

(8) 
(9) 



As standard in perturbation theory we assume that the eigenvalues and 
eigenfunctions of eq. Q can be expressed as power series of the perturbation 
density 5p, i.e. that 



k=0 



and 



<^ n (x) = J2^(x)r 1 k 



(10) 



(11) 



k=0 



Clearly E { ° ] and $l 0) (x) are the energies and the solutions for a uniform 
string of density po; we also call e n the energies of a uniform string of unit 
density and use the Dirac bra-ket notation to represent the modes. 

Readers familiar with perturbation theory for stationary states in Quantum 
Mechanics, i.e. with the Rayleigh-Schrodinger perturbation theory (RSPT), 
may easily spot the similarities and the differences between the typical imple- 
mentation of RSPT and the present case: the approach that we follow here is 
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the standard one used in the treatment of time independent perturbations in 
quantum mechanical problems, although in a typical quantum mechanical prob- 
lem the Hamiltonian is polynomial in the coupling (most likely a linear function 
of g). In the present case the perturbation is the density 5p(x) and it appears 
in a non-polynomial form in the symmetrized operator of eq.Q. 

As a result, at a given order in 77 there will be contributions stemming from 
the operators of different orders O n allowed to contribute to that given order 
and therefore there will be a proliferation of terms as higher orders in r\ are 
considered. 

The first few corrections to the energies of the string have been worked out 
explicitly in ref. [l], 0] and read: 

Po4 0) = e„ (12) 



po 
Po 



4 1 ) = -e n (n\6p\n) (13) 
p(2) / ix 1 \2 , 2 \ - (n\Sp\k) 2 

E i> = e n (n\Sp\n) + e„ ^ ~ ( 14 ) 

k^n £ ™ 6 ' 

= -e n (n\6p\n) 3 + e 3 n (n\Sp\n 



(n\5p\k)' 



q 2/ ir 1 \ V- ( n \Sp\k) 
3e n (n\6p\n) ^ 



k^n 

2 



U) n k 

k^n 

(n\5p\k)(k\5p\m)(m\5p\n) 

'"2^,2^ ' ( 15 ) 



where e„ = n 2 n 2 /p and iu nk = e n -e k . 

Here the reader will easily recognize that the direct terms at each order, i.e. 
the terms involving matrix elements of Sp between the same states, appear to 
originate from a geometric series. These terms may be resummed into 

(n\p(x)\n) 

which for n — > 00 would lead to an asymptotic behaviour 

lim E n \^ dlrect ^ a —j^ . (17) 

n ^°° J_£ p{x)dx 

In two dimensions, when p(x) is the density obtained from a conformal map 
which transforms an arbitrary domain into a square, the denominator of this 
equation provides the area of original domain, thus leading to a alternative 
derivation of Weyl's law, based on perturbation theory [5j. The validity of 
this result however requires that the direct terms dominate asymptotically the 
terms in the perturbation series which mix the different states. A specific case 
discussed in ref. [H concerning the small deformations of a square, has been 
proved to fulfill this requirement to second order. 
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In one dimension the situation is different: for this case we have proved in 
ref. [l[ that for n — > oo 



IX) 



(18) 



This equation is consistent with the asymptotic behaviour which is obtained 
using the WKB method 0], 



(19) 



where (^/p) = ^ J^ L y / p(xjdx. 

One could regard this equation as Weyl's law for an inhomogeneous string. 

We will now apply the variational principle to obtain useful upper bounds 
to the energy of the fundamental mode of a string with arbitrary density. First 
we observe that the spectrum of the hermitian operator O corresponding to 
this string is positive and bounded from below; if we let ^(x) be an arbitrary 
function obeying Dirichlet boundary conditions at the ends of the string then 
the expectation value of O in this state is an upper bound to the exact energy 
of the fundamental mode: 



(20) 



The proof of this statement is straightforward once ^f(x) is expanded in the 
basis formed by the eigenstates of O. 

Working to lowest order we use the fundamental mode of a uniform string 
of density po as a variational ansatz 



1 

— sin 



(x + L) 



2L 



(21) 



and obtain the bound 

(Oh = 



+L n 2 CQS 2 



-L ^p{ X) 

+L cos 2 (ff) p"{x) 
_ L 2Lp(xf 



dx + 
dx 



+ L 



™(!!)cos(fg)p'(*) 



L 3cos 2 (ff)p'(x) 2 



dx 



-L 



ALp{xf 



dx > Ei (22) 



This bound is expected to be useful only for slightly inhomogeneous strings. 



3. WKB perturbation theory 

We now describe a different approach to the problem of a string of variable 
density. This approach has been first used by Lindblom and Robiscoe in ref. @ , 
where the inhomogeneous Helmholtz equation was cast in the form of an equiv- 
alent Schrodinger equation on the finite domain and approximate expressions 
for the eigenvalues of eq. ([T]) were obtained. 
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The path that we follow here is similar to that of ref. 6], although our 
primary interest lies in the development of a perturbative expansion for the 
eigenvalues and eigenfunctions of eq. (HJ. 

We first define the functions 



a(L) 



p(x) 



1/4 



sin 



a(x) 



a(L) 



n = l,2,... 



(23) 



where 



a{x) 



\fpiy)dy 



(24) 



After making the change of variable u = a(x)/a(L) it is straightforward to 
see that that these functions are orthonormal on [— L, L], i.e. that 



+ L 



-L 



4>n{x)4>m{x)dx = S n 



(25) 



This should not be a surprise since the change of variable relates these functions 
to the eigenfunctions of a uniform string. Moreover, the <j) n {x) are closely related 
to the WKB functions (see for example eq.(10.1.31) of Bender and Orszag, 
ref. at pag.490). 

Let us now consider the action of the symmetrized differential operator O 
which has been introduced in the previous section on any of the functions given 
above; with elementary algebra we see that 



d(j> n (x) 



W Ap{x)p"{x) - bp'jxf 



lQp{x) 



(j>n{x) 



(26) 



For the highly excited states (n ^> 1) the second term in the parenthesis 
is neglegible with respect to the first and the (j) n {x) become increasingly good 
approximations to the exact eigenfunctions of O, as usual in the WKB approx- 
imation. 

The <f> n are eigenfunctions of the hermitian operator Q: 



b n (x) = 



O- 



4p(x)p"(x) - 5p'(x) 2 



16p(x) 



<t>n(x) 



n 2 n 2 



(t>n{x) , 



(27) 



whose eigenvalues coincide with the leading WKB expression for the eigenvalues 
of eq. (Trj. 

In special cases the <j)n{x) provide the exact solutions to the full problem: 
this happens when the density is such that 



4p(x)p"{x) - bp'jx) 
16p(x) 3 



K , 



(28) 



where k is a constant. 
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The general solution to this ordinary differential equation has the form 

256c 2 

P[X) ~ (cj (c 2 +x) 2 - 256 K ) 2 ' ( ' 

where c\ and c 2 are integration constants; the density discussed by Borg in [7] 
is a special case corresponding to 

16a 2 1 

ci = , c 2 = L + - , k = . (30) 

1 + a a 

With the purpose of implementing a perturbative expansion we define the 
operator 

O v = Q + r 1 ^pte) 3 =Q + r)V(a(x)) , (31) 

where r\ once again is a non-physical parameter used for power counting (0\ = 
O) and V(a(x)) is a potential defined in terms of the density. We assume that 
the eigenvalues and eigenfunctions of O can be expressed as power series in 
77, as done before in eqs. (fT0|) and (fTTj) for the DPT. Unlike for the DPT we 
notice that here the density p(x) is not bound to be a perturbation and that O v 
is an operator which depends linearly on 77. Because of this last property the 
corrections for the energies of the string in this approach have a simpler form 
than in those in the DPT. 

We write explicitly the first few terms: 

= 02) 
= (n\V\n), (33) 

(34) 



(2) a_{Ll_ v \(n\V\k)\ 

n " TT 2 ^ (n 2 - k 2 



^(^) 4 V- V- (n\V\k)(k\V\l)(l\V\n) 



in 2 - k 2 )(n 2 - I 2 ) 



Using eq. (I3T1) we may cast these expressions directly in terms of the operator 
O, taking into account that (n\0\k) = (n\V\k) for k 7^ n. Notice that the second 
order correction to the energy of the fundamental mode is always negative as in 
the standard RSPT. 

Similarly for the eigenfunctions we may write 

4>W{x) = <j> n {x) (36) 



7 



a(L) 2 ^ (k\V\n) 
n 2 — k 2 



E 



<t>k{x) 



EE 



(k\V\j){j\V\n) 
(n 2 — k 2 )(n 2 — j 2 ) 



<i>k(x) 



G ^ A I IT/I >r ( k \ v \ n ) , , n 
— —(n\V\n) ^ rra ^H 

_ (n |T/| n) ^ ^ (n2 _ fc2 )2 

<x(£) 6 , .... , 2r (&|V» , 



a(L) 6 



EEE 



{n\V\l){l\V\n){k\V\j) 
(n 2 — k 2 )(n 2 — j 2 )(n 2 — I 2 ) 



4>k(x) 



k^n j^n l^n 



EE 



(n\V\l)(3\V\n) 2 
(n 2 — k 2 ) 2 (n 2 — j 2 ) 



■4>k(x) 



(37) 



(38) 



(39) 



The expressions that we have written above are equivalent to those found in 
ref. @; in particular eq.(2) of that paper is just the energy of a string calculated 
to first order in this scheme. 

We will now derive asymptotic expressions for the energies of a string of 
variable density using these perturbative formulas and compare them with the 
asymptotic formulas obtained by Lindblom and Robiscoe. 

As we have already mentioned, eq. ([32]) is just the leading term for the energy 
of an inhomogeneous string obtained with the standard WKB method [H . The 
differences with the WKB results appear at first order, in En] with a simple 
change of variable we may cast eq. ([33| as 



M 1 ) = 



<L) Jo 



sin 



a(L) 



V{a{x))da . 



(40) 



Eq.(2) of ref. [6| contains this expression. 

We may now use the results of |Appendix A to write E,P in the equivalent 
form 



oo / . \ 2k-\-l 

4 1} = (V) - E ( g 2 ; n)2fc+2 (- 1 ) fc \v {2k+1 \«{L)) - T/( 2fe + 1 ) (o) 



fc=0 



(41) 



where (V) = ^ f* L) V(a(x))da and V^ 2k+l \a) = ^^£11. Notice that 
no approximation is made in deriving this expression. 
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The first term in this expression is independent of n and clearly dominates 
for ri 3> 1. We may write it explicitly as 



(V) 



1 



a(L) 



+ L 



4p(x)p"(x) - 5p'(xf 
16p(x) 5 / 2 



dx 



or equivalently as 



(V) 



£tp{xy/Hx 



(42) 



(43) 



This term agrees with the first correction to the energy that can be ob- 
tained within the WKB approximation (in the WKB approximation however 
the " quantization" of energies follows from enforcing the boundary conditions 
at one end of the string) and corresponds to the coefficient of eq.(25) in 
ref.Q. 

We now consider the contributions of order 1/n 2 , which originate both from 
E$P and We write eq.flSJ to this order as 



= (V) 



V {1) (a{L))-V {1) {Q) 



+ o 



(44) 



where 



= dV(a(x)) dV(a(x)) 1 

do dx ^/pjx) 

p"'{x) 15p'(x) 3 9p'(x)p"(x) 



Ap(x) 5 / 2 16p(x) 9 / 2 8p(x) 7 / 2 



(45) 



A contribution of order 1/n 2 stems also from . Using the expression for 
the off-diagonal matrix element (n| V|fc) obtained in | Appendix A we may write 



a 2 (L) ^ {n\V\kf 



E 



-k 2 



77/ 



ky^n 

2 / T \ oo oo 



a 2 (L) 



fe^n 3=0 (=0 
1 



(L) 



1 



1 



7r 2(j+/)+4 n 2 _ k 2 
1 



[k - n) 2 J+ 2 (k + nfi+ 2 ) \(k- n) 2l + 2 (k + n) 2l + 2 
{-l) k + n V {23+1) {cr{L)) - V W+1) {Q) 

(-l) k+n V^ +l \a{L))-V {2l+l \Q) 



where no approximation has been made. 
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This expression may be cast in the more compact form 



M 2) 



^2 Z^Z^ ^ 7r 2(j+i)+4 
j=0 1=0 

{s$\n) [v^ +1 \a(L))V^(a(L))+V^ +1 \0)V^ l+1 \0) 



- S 



(ii) 



(n) [v^ +1 \a(L))V (2l+1 \0) + V^ +1 \0)V {2l+1 \a{L))\ j(4p) 



where we have defined the series 



Y — 



1 



1 



n 2 -k 2 \{k-n) 2 i+ 2 (fc + n) 



2j+2 



1 



1 



2;+2 



(k-n) 

(_!)fc+r. 
Z^ „ 2 _ 



(k + n) 



21+2 



n 2 -k 2 \(k-n) 2 i+ 2 (k + n) 2 J+ 2 



(k - n) 2l + 2 



(k + n) 2l + 2 



We are interested in calculating the dominant contributions in these expres- 
sions for n — »■ 00: 



+ 



c(2(2+j + q) „ r 1 

2n 2 



C (2(2+i + 0) r I, 2 -2(2+i+0 
~2 + 





' 1 " 


+ 


n 3 



C(2(2+j + i)) 



2n 2 



1 OO 

n 2 f-^ 



(2fc) 2 ( 2 +J+0 



O 



(47) 



(48) 



Therefore we may write: 



E ( 2 ) 



21 t\ °° °° 2(j+i)+2/r\ 

n 2 n 2 Z^Z^ ^ n 2(j+l)+A 
3=0 1=0 



1 1 



V^ j+1 \a(L))V^ l+1 \a(L)) + (O)^ 2 '^ (0) 



(2 ^ (jfe)2(2+j+i) 

+ T/^'+ 1 )( ( 7( J L))T/( 2i + 1 )(0) + V^ +1 \Q)V {2l+1 \a{L)) 
1 



£ 

o 



^ (2fc) 2 C 2 +-7+0 
1 



v«»' +1 >(<r(L))v< 2,+1 >(0) + v^ +1 Ho)v^ 2l+1 \a{L)) 



(49) 
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Let us define the series 



3=0 



which is related to the series 
G(v,<t ) = 



j2j+l 



3=0 



(2j + l)! 



u2j + l y (2 J + l) ((To) 



by a Borcl transform 



F{v,a a ) 



[V(a + iv) - V(a a - iv)] 



(2) 

We can therefore cast E n directly in terms of T as 



a(L) 



2 00 



T 



k=l 

2nk ' 



T 



2k 2 
a(L) ) T 



irk 

<r(L) 
2nk 



a(L) 



o(L) 
irk 



,0 



o 



To the best of our knowledge this result was not known. 

Therefore, the exact expression for the energy to order 1/n 2 is 



a(L) 2 
a(L) 2 



n 2 + (V) + 4j 



oo ^ 

^ 2k 2 



I 



k=l 

a{L) 



4tt 2 
nk 



,a(L) 



V (1) (a(L))-V (1) (0) 
a(L) 



2nk 



T 



O 



irk 



(50) 



(51) 



(52) 



(53) 



In the next sections we will explicitly check this formula working on specific 
examples. 

To make contact with the results of Ref. [6[ we consider a simple approxi- 

(2) _____ 

mation to E n , which is obtained restricting the sums over j and I in eq. (|46l) 
to j = I = 1_| in this case we have 



180n 2 



O 



(54) 



x We will soon discuss examples where this series is divergent. 

2 Clearly this formula is not expected to be accurate, unless the contributions to the term 
which goes like n~ 2 stemming from -E^ 2 ' are particularly small compared to the one stemming 
from s£ . 
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7 



1440n 2 

and the energy is approximately given by 

E (2) 



+ o 



(55) 



180 



(vW(a(L))) + (V«(0))' 



(56) 



To allow the comparison with the eq.(15) and (20) of [6j we write the ap- 
proximate analytical expression for the energies of the inhomogeneous string 
using eq. (|56"|) (i.e. restricting the double series in the expression for the 1/n 2 

(2) 

coefficient inside E„ to the first term): 



+ 



cr(L) 2 n 2 [ 47T 2 

ct^L) 4 rFW(a(i)) 2 + FW(0) 2 



180 



7W(ff(L))-VW(0) 



^V^(a(L))V^(0) 



o 



1 



(57) 



There are some important considerations to make: to start with, it should 
be easy to convince oneself that the asymptotic expansion of the energy is an 
even function of n and that contributions to a given order 1 / n 2p can only come 
from the first p perturbative orders. This happens because the non-diagonal 
matrix elements go to zero as n — > oo and because the WKBPT contributions 
to a given order p contain p of these matrix elements (an example to this is 
given by the calculation that we have just performed, where the corrections 
to the energy decaying as 1/n 2 only come from and E^). This means 
that the expression that we have written above is exact to order 1 /n 2 and that 
the evaluation of higher order terms can be done by using a larger (but finite) 
number of perturbative terms. 

The second point concerns the relation of our expression with the one ob- 
tained by Lindblom and Robiscoe, who estimate the second order contributions 
using the mean value theorem for derivatives. The two expressions are not 
equivalent since the results of ref. [6[ are obtained assuming that the density 
is a slowly varying function as stated in their eqns.(ll) and (12); we have not 
made this hypothesis on the density in our calculation. Morevover the term 
V{x) in the equation (I26[) is always a perturbation for n sufficiently large pro- 
vided that V(x) is not singular on the domain, as previously mentioned. It is 
straightforward to see that our result of eq. (|56"|) reduces to the result of ref. Q 
in the case where V^- 1 ' is a symmetric function of x: in this case we have that 
VW(ct(L)) = yW(0) and 2/180 + 7/720 reduces to the factor 1/48 of ref. @. 

We will now use the WKBPT to obtain bounds on the energy of the fun- 
damental mode of a string of arbitrary density, using the solutions in eqs. (|36p 
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and (f3"T|) to build a suitable variational ansatz. As we have seen in the previ- 
ous section, for an arbitrary (normalized) ansatz ^(x), the Rayleigh quotient 
provides an upper bound to the true energy: 

(6)* > Eh. • 

In this case we will use the wave function for the fundamental mode obtained 
using WKBPT to a finite order as a variational ansatz; the simplest approach 
corresponds to choosing the lowest order expression as the ansatz 



M>(z) = ^°)(z) 
which allows one to obtain the bound 
r 2 2 



(58) 



(0)» = 



a{Lf a{L) 



+ L 



sin 



it<j(x) 

Mlj 



V(x)\/p(x)dx 



(59) 



which is just eq.(2) of ref. [6j for n = 1, although the authors of that paper do 
not mention this property explicitly. 

The first term in this expression is the leading WKB expression, which typi- 
cally is not very accurate for the fundamental mode, and it is always positive; the 
second term, on the other hand, depends on V{x) and it does not have a definite 
sign. In the WKB method this last term reduces to f_ L V(x) \J p(x)dx, 

t{x)~ 



having substituted the sin ~njy with its average value on the interval. Un- 
fortunately, the WKB expression neither provides a bound nor it is accurate, as 
we will later see in a specific example. 

If we like to obtain a stricter upper bound for the energy of the fundamental 
mode of our inhomogeneous string we may use the ansatz (not normalized) 



*( X ) = ( j ) { 1 o) (x)+ ( i ) r>(x) , 



(60) 



which corresponds to the solution for the fundamental mode calculated with 
WKBPT to first order, restricting the sum in the first order contribution to the 
first N — 1 terms: 



In this case 



where 



(0>* 



<r(L) = 



a(L) 



2 ^(k\V\l) 



ML) 



^ 1 

k=2 



<t>k{x) 



6i 1)2 



4 N 



> E x 



Y- 

^ (1 

k=2 



\\V\k) 2 
(1-fc 2 ) 2 



(61) 



(62) 



(63) 
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a(L) 



a(LY 



2 N 



^ (i-fc 2 ) 2 ' 



fc=2 



<r(L) 



2 N 



Y- 



(1-fc 2 ) ' 

fc=2 V ' 

{\\V\k)(\\V\l)^l\V\k) 
(l-/c 2 )(l-Z 2 ) 



(64) 



(65) 



We will also apply this second bound later to specific examples. 



4. Improved WKB perturbation theory 



We will now describe an alternative implementation of WKBPT which is 
suitable for describing specific parts of the spectrum. 

Let us consider an unphysical density p(x): by unphysical density we mean 
a positive function defined on the domain of the string (p(x) > for \x\ < L) 
and different from the physical density of the string, p(x) ^ p(x). 

We define (j) n {x) the functions corresponding to eqn. f[23]) for p{x) — > p{x) 



4>n(x) 



a(L) 



^(x) 1 / 4 sin 



a(x) 



a(L) 



n = l,2,. 



(66) 



where 



(67) 



The <j) n {x) are eigenfunctions of an operator Q, obtained form Q with the 
substitution p{x) — > p{x) and they are isospectral with the 4> n {x) if a(L) = a(L). 
If we apply O to 4> n (x) we see that 



0(f) n (x) 



n 2 7r 2 p(x) 



n 



4>n(x) 



(68) 



i.e. that only the functions corresponding to p(x) = p{x) are eigensolutions 
of O for n — > oo. On the other hand, one can pick a specific density p(x) so 

that a given state of Q is also eigenstatc of O. We can think of this density 
as an effective density, which will reduce to the physical density for states with 
n — > oo. 

This density may be obtained as follows: suppose that ^ n (x) is the n th 

eigensolution of O, and let 4> n (x) be the n th eigensolution of Q. We want to 
find the density p(x) for which 



^ n (x) = 4> n (x) , 



(69) 
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which is equivalent to the equation 



a{L) 





2n-n<j{x) 


sin 


a(L) 



2nir 



(70) 



Once cr(x) is known, the density is obtained using the relation a'(x) = \/ p(x). 

The procedure outlined above explains how the effective density for a given 
state of O can be obtained once that the eigenfunction is known explicitly; in 
a typical situation however the eigenfunctions of O are not known and need to 
be calculated. In this case one must resort to an alternative procedure. 

If we are interested in describing the lowest mode of a string with density 
p{x), we may use the variational principle to obtain the effective density p(x) for 

this mode: in this case the eigenfunction corresponding to the lowest mode of Q 
can be used as a variational ansatz and therefore one can obtain the inequality 



p{x 



,1/4 



°( L ) J-L ^/pjxj 



TT&(x) 



dx 2 



p{x) 



1/4 



1Ta(x) 



dx>Ex, (71) 



The LHS of this inequality is just the expectation value of O in the funda- 
mental mode of a string with density p(x) □• If p(x) depends on one or more 
parameters which can be varied without p(x) changing sign on all domain, then 
one obtains an optimal approximation to the fundamental mode of O in corre- 
spondence of an absolute minimum of the LHS of the inequality. The density 
obtained using this values of the parameters approximates the effective density 
for the fundamental mode of O. 

The bound discussed above can still be improved: Gottlieb has shown in 
ref. Q that a string of length 2L centered in the origin and with density 



p(x) = £'(*) 2 p(£(x)) 
where a > —1/2L is a real parameter and 

1 + 2aL , 
1 + a(x + L) 



(72) 



(73) 



is isospectral to a string of the same length also centered in the origin and with 
density p{x) 0. Using this property one may express the bound (|7"Tj) for the 

operator O = * ( — -j-^ ) * . By choosing a suitably one may obtain a 

ypO) V ax I y p{x) 

more accurate bound for the energy of the fundamental mode of the string. 

In the case of arbitrary excited modes of O the approach to obtain the 
effective density for these modes is different. Let us consider an arbitrary trial 



3 Notice that (O) is invariant under the redefinition of the density p(x) — > cp(x), where c 
is a constant. 

4 We may easily prove this property by proving that V is invariant under this transforma- 
tion. 
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function & n (x) — (x\^ n ), with \& n (±L) = and with n—1 nodes in the interval 
x € (-L, L). We wish to obtain accurate approximations to the eigenfunctions 

of the hermitian operator O — } \—4-*\ , 1 For this reason we define 

the quantity 

En = <*„|0 2 |*„) - (*n|0|*„) 2 > , (74) 

where the equality holds only when is a solution of O. 

Provided that we vary *Sf n (jx) without changing the number of nodes in 
(-L, +L), the optimal approximation to the n th solution of the string falling in 
this class of functions will be the one which minimizes £„. We can achieve this 
goal using a trial wave function corresponding to the n th solution of a string of 
arbitrary density p~(x), which depends on one or more arbitrary parameters: 



= ~ K{x) = fa?/* S in . (75) 

By varying the parameters in p(x), with p(x) > on the interval, we modify 
the trial solution, without changing the number of nodes and therefore the 
minimization of the corresponding S n will yield the optimal approximation to 
the mode that we want to calculate. The corresponding p{x) will approximate 
the effective density for this mode. 

From a computational point of view the calculation of E n as given in eq. (|74j) 
is not efficient since the difference of small quantities may lead to large numerical 
errors and even to negative values of £„. A more efficient calculation can be 
done using the completeness of the orthonormal basis to which <j> n {x) belongs 
and writing 

oo 

E n = J2(n\0\k) 2 , (76) 

where (x\k) = (f>}.(x). In this case £„ is the sum of positive quantities and 
therefore the numerical errors do not affect the sign of E n . 

Provided that the matrix elements (n\0\k) fall off sufficiently rapidly away 
from the principal diagonal, the series may be truncated to a finite number of 
terms, without introducing large errors: 

4 n) « E HO|fc) 2 , (77) 

kmin ,k^n 

where k m i n = max(fc — m, 1) and k max = k + m. 

We are now ready to describe a new perturbation scheme which uses the 
effective density of a given state of O and the set of functions related to this 
density, <p n (x), as the "unperturbed" basis. We write the operator 

Or, = Q + v \0 - Q\ = Q + rjW , (78) 
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where Q is the operator obtained from the hermitian operator corresponding to 
a fictitious string of density p(x), 



O 



1 



1 



by means of the decomposition 

d = Q + V(x) , 

with 

4p(x)p"(x) - 5p'(x) 2 



V(x) 



16p(x) 



(79) 



(80) 



(81) 



Notice that when p{x) is chosen as a constant density po we have that V = 
and the operator Q = O. In this case one may write 



Or, 



po dx 2 



1 



1 





y/p{x) dx 2 y/p{x) Po dx 2 



(82) 



The Density Perturbation Theory (DPT) discussed early is recovered after 
expressing p(x) — po(l + Sp(x)), with po the average value of p(x) over the 
interval, and expanding in powers of 5p(x): therefore DPT may be regarded as 
a special case of the present scheme. 

Going back to the general expression, we point out that 77 is a parameter 
allowing power-counting and that O^—i = O; we treat W as a perturbation 
and obtain the perturbative expressions for the energies and eigenfunctions of 
a string of density p(x) as[3 



n 

4 2) = 
M 3) = 



a(L) 2 ' 
(n\W\n) 



a(L) 2 ^ \ (n\W\k)\ 2 
(n 2 - k 2 ) 



E 



(83) 
(84) 
(85) 



HL) 4 



HL) 4 



EE 

k^n l^n 



(n\W\k)(k\W\l)(l\W\n) 



i\W\n) ]T 



(n 2 - k 2 )(n 2 - I 2 ) 
n\W\k) 2 



(86) 



k^n 



5 Here we call |n) the eigenstates of Q, thus implying that (x\n) = (j> n (x). 
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and 



<t>£H*) 



4>n(x) 

a(L) 2 ^ (k\W\n) ~ 



HL) 4 



a(LY 



7T 



E 

k^n 



EE 

k^n j^n 



i 2 -k 2 



(k\W\j)(j\W\n) 
(n 2 — k 2 ){n 2 — j 2 ) 



4>k(x) 



(n\W\n) 



k^n 



(k\W\n) 7 , x 
{ n 2 -k 2 ) 2Mx) 



7T U 



+ -^-H^l") 2 E ct ^W 

k^n 

c^(L) 6 vvv H^|Q(/|^|n)(fc|^|j) - 

^6 Z^Z^Z^( n 2_ fc 2)(„2__ 7 -2U n 2_m9'fcW 

? ( L ) 6 / IWI \ \^ (n\W\l)(l\W\n) - 



k^n l^n 



(n 2 - k 2 ) 2 (n 2 — I 2 ) 



HL) 6 ^ ^ (n\W\l)(j\W\n) 2 - 



k^n j^n 



(n 2 — k 2 ) 2 (n 2 — j 2 ) 



(87) 
(88) 



(89) 



(90) 



Having chosen p(x) to approximate the "effective" density for the n th excited 
state of O we expect that this perturbation scheme will be helpful to obtain 
precise approximations to this state and to states which fall near since in this 
basis the n th state is almost decoupled from the remaining states. 

We will now discuss the asymptotic behaviour of the energies calculated to 
first order: 



£i 0) 



2 



a{Lf 
+L p(x) lfi . 



n\W\n) 



nira(x) 



a(L) 



dx 2 



1/4 



&( L ) J-L y/pjxj 

For n — > oo the leading contribution goes as 

~a(L) 2 f +L p{xf/ 2 



mra(x) 



dx . 



E<P+E<P 



a{Lf 



a{Lf 



p{x) 



dx 



(91) 



Deviations of the term inside the parenthesis from 1 imply violations of the 
correct leading asymptotic behaviour. 
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5. Iterative methods 



Another approach which can be used to obtain precise approximations to 
the solutions of the equation for an inhomogeneous string has been recently 
discussed in ref. In that paper we have proved three theorems which allow 
one to build iteratively the solutions of the fundamental and excited modes 
of the string starting from a trial solution, fulfilling the appropriate boundary 
conditions. 

In particular, the main theorem, Theorem 1, may be view as a generalization 
of the well-known power method of linear algebra to the inverse operator of 
the inhomogeneous Helmholtz equation: the sequence of iterations produces 
approximate solutions which provide a sequence of increasingly tighter bounds 
for the energy of the fundamental mode when the expectation value of the energy 
is calculated. The only limitation in the application of these theorems is given 
by the ability of performing the relevant integrations explicitly. 

We state now the three theorems, which are proved in ref. [1]. Applications 
of these theorems will be considered in the following section. 

Theorem 1. Let t;o(x) be an arbitrary function defined on the interval (— L, L) 
(being S,q(x) arbitrary we assume that it has nonzero overlap with the true fun- 
damental solution). For n — >■ oo the sequence of functions 



£ n (x) ee \/p(x) I ^dy 



dzy/p(z)£ n -i(z) 



(92) 



with K n = -jj- J_ L dy J_ L dz^J p(#)£ n _i(z), converges to the fundamental mode 
of the Helmholtz equation (QJ) with Dirichlet boundary conditions, £ n (±L) = 0: 

* o (a0 = lim Cn(x) , (93) 

where (x) ee . 



Using this theorem we may write the bounds 

+Z — • • • — +1 — 1 ' ) 

Lh £l( X ) dx I-L Zn(x)dx 

where E\ is the exact energy of the fundamental mode of eq. (TTJ) . 

Theorem 2. Let ee j^o 1 ' (a:), . . . j^o^C 2 -)} ^ e a se ^ °f arbitrary functions 
defined on the interval (-L, L). We assume that these functions have a nonzero 
overlap with the first N solutions of the inhomogeneous Helmholtz equation. For 
convenience we assume that the functions are ordered as 

Q j^\x)d^\x)dx < < J I : L ^ N \x)6^ N \x)dx 
j L _ L &\x)*dx £ L ^\x)Hx 
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Consider the sequence of functions 



= Jpjx) / dy - / dz^pjz)$\z) 



(96) 



with 



Let 



Jl S-l d V I-l dz Vpiz)£,u ] (z)andj = l,...,N 



U), 



(97) 



and = ^[^(x) , . . . , ^[ N \x)^ . We call EE^™-* the set of functions obtained 

after n iterations. Then, for n — > oo, 3^ n ' converges to the N lowest eigenfunc- 
tions of eq. J7J). 

Theorem 3. Let T]o(x) be an arbitrary function fulfilling Dirichlet boundary 
conditions on [— L, L'], with —L < L' < L. The sequence of functions 



r) n (x) = \J p{x) I dy 



dzy/p(z)r] n -i(z) 



(98) 



with k n — L ^ L , f_2 dy J^ L dzyj p(z)r) n -i(z) , converges to a solution of equa- 
tion 0) for n — > oo if L' is chosen so that rj n (L) = 0. 

This theorem can be used to obtain estimates for the energies of the excited 
solutions of eq. ([!]): 



E k 



f- L m {x)r} { n\{x)dx 



(99) 



dx 



where Ef. is the exact energy of k th state and r/n (x) is the function obtained 
after n iterations and corresponding to the value of L' obtained from the k th 
solution to the equation rj n (L) = 0. Notice that the Ek for k > do not provide 
upper bounds to the exact energies. 



6. Some applications 

In this section we consider some applications of the approaches described in 
the previous sections. 

6.1. The ratio of the first two eigenvalues of an inhomogeneous string 

A first application concerns the calculation of the ratio of the first two eigen- 
values of a string of variable density. Huang has proved in ref. @ that for strings 
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with concave densities the ratio EijE\ is minimized when the density is con- 
stant, whereas for strings with convex densities the ratio is maximized when the 
density is constant. 

We may check what DPT and WKBPT can tell us concerning this problem. 
Working to first order and assuming p(x) = 1 + 6p(x) with f^ L 5p(x)dx = 0, 
we obtain the approximation: 

J«4[l-<2|«5p|2> + <l|Ml)] , (100) 

which holds for \Sp\ <C 1. 

To recover the result of Huang we rewrite our formula as 

g«4[l-<2|p|2> + <l|p|l>] . (101) 

When the density takes larger values at the ends of the string one has that 
(2|p|2) — (l|p|l) > 0, since the solution for the second mode takes larger values 
(in absolute value) in this region with respect to the solution of the first mode: 
in this case the ratio of the first two eigenvalues is maximized when the density 
is constant; conversely, when the density takes larger values at the center of the 
string, then (2|p|2) — (l|p|l) < and the opposite behaviour takes place. 

Let us now consider the same problem using WKBPT: we work to first order 
and obtain the ratio 

| 2 -«4+^ ({2\V\2)-4{1\V\1)) . (102) 

E\ TX Z 

As before we assume p(x) = 1 + 8p(x), with 8p(x)dx = 0: assuming 
that \5p\ C 1 we may approximate V(x)o~'(x) ~ \8p"(x) and obtain a simple 
expression for the ratio: 



E X 7T 2 



+l r 

1 : , \ 1 —1— 1 . \ \ 1 i 1 i ■ / 1 \ \ f j 

Sp (x)dx 



4 sir —— — - - sm^ ' 



2L 



The term in parenthesis is positive and peaked at the center of the string, 
therefore the sign of the correction to the ratio of the energies depends essentially 
on the sign of 8p"(x) in proximity of x = 0, confirming the analysis done using 
the DPT. 

6.2. An exactly solvable string 

Horgan and Chan give in ref. flo| explicit examples of strings with nonuni- 
form density whose frequencies and solutions can be calculated analytically. We 
will use one of these examples to test the accuracy of the asymptotic formulas 
for the energies obtained before. 

The density of this string is 

P{X) = 2(a(a+ ( 2)(2x + l) + 2) ' (1 ° 3) 
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Table 1: Coefficients of the fit of the energies of the string with density (1 1 03 I t with a = 1. The 
third column contains extrapolations obtained from the values of the fit. 





Fit 


Extrapolation 


|Fit — Extrapolation! 


A 2 
A 3 
A± 
A 5 


9.8696044010893586188 
0.3750000000000000000 
-0.0326523345722377584 
0.0091705849049749141 
-0.0058099597602144426 


it' 2 

3 

165 

«443 
14fiSnnfi47r6 


1.5 x 10" iU2 
1.1 x 10~ 93 
3.5 x 10~ 85 
6.8 x 10~ 77 
9.1 x 10~ 69 



where a > —1 is a real parameter. The string has unit length and it is centered 
in the origin. 

Horgan and Chan have shown that the frequencies of this string are obtained 
from the zeros of the trascendental equation 

Ji(A)Yi((a + 1)A) - Yi(A)Ji((o + 1)A) = (104) 

as 

uj = aX. (105) 

Here J and Y are Bessel functions of first and second kind respectively. 

It is therefore possible to calculate a large number of energies (E = oj 2 ) for 
this string with high accuracy, and therefore estimate the asymptotic behavior of 
the energy by means of a simple fit with the known leading asymptotic behavior 
of the energy 

E n « Am 2 + A 2 + A a /n 2 + ... (106) 

Wc have chosen a = 1 and we have calculated numerically the first 10000 
energies of this string with high precision (each zero has been obtained to 200 
digits). We have then used the energies from 5000 to 10000 to obtain a fit as in 
eq. (|106|) involving the coefficients A\, . . . , An. The first 5 coefficients obtained 
with this fit are reported in the second column of Table [TJ 

In the third column of this table we report the explicit expressions to which 
the numerical coefficients in the second column seem to converge: the calculation 
of the first two coefficients A\ and A 2 using WKBPT is straightforward and 
confirms the corresponding results in the third column. 

A more interesting test is provided by the calculation of A3 using WKBPT; 
for this problem we have that <j(L) = 1 and V(cr) = ^jj^p ■ As we have seen 
the coefficient A3 receives contributions both from the first order and the second 
order term in the perturbative expansion: 

A 3 = 4 X) + 4 2) - (107) 
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where A^ (A^) is the contribution from first (second) order WKBPT. 



The evaluation of the first term is trivial: 



.4 



(i) 



At: 2 



V {1) (a(L))-V (1 \0) 
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-0.03324601338 . (108) 



Notice that this term accounts for most of the value of A 3 reported in the table. 
Let us now consider the second contribution; in this case 



Q(v, ct ) = -- [V{a + iv) - V(a - iv)} 
3(<j + l)v 



2((<T + i) 2 + v 2 y 



(109) 



where in this case Co can take the values and a(L) = 1. 

The calculation of T can be performed analytically and one obtains H: 



T(v,a Q ) 



[2fi Ci(fi) sin(fi) - 2fl Si(ft) cos(ft) + ntt cos(ft) -fl2]0) 



8£lv 2 

where S7 = +a ° and Si and Ci are the sine and cosine integrals. 



Having obtained an explicit expression for T we can then evaluate the A 3 
coefficient, which has been cast in the form of single series as: 



(2) 



.4 



(2) 



a(L) 2 



^ 2k 2 



k=l 



nk 



1 2 



2irk 



fe=l 



where 



(111) 



Ik 



1287T' 



{U{-l) k k{-K - 2Si(fe?r)) - 27rA:Si(2fc7r) + n 2 k - 3) 2 
-27r/fcSi(2fc7r) + 7r 2 A:-l) (-47r/fcSi(4A:7r) + 27r 2 fc-l)) . (112) 



The series for A 3 may be approximated by a partial sum over the first iV 
terms in the sum: for instance when we restrict the sum to the first N — 5000 
terms we have 



5000 



A 3 w A { 3 > + Ik ~ -0.032652334572241609068 , 



k=l 



(113) 



6 Notice that 



{ 0) 2 (l + CT0 )2.+3 



and therefore the series defining J^(v, (Tq) is divergent. 
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which agrees with the result in the table to 14 digits. 

A much better estimate may be obtained taking into account the asymptotic 
behavior of the 7^ for k —> 00: 



711 



Ik ~ 

and therefore we may write 

l(2) ^ 79 



5127T 6 fc 4 



4 2) =E^ = 



k=l 



5120tt 2 



E 

fc=l 



711 



7fc 



5127T 6 fc 4 



(114) 



(115) 



If we limit the sum to the first 5000 terms as before we obtain a more precise 
estimate for A3: 



A, 



.4 



(i) 



79 



5120tt 2 



5000 

E 

k=l 



Ik 



711 



5127T 6 fc 4 



-0.032652334572237758375 , 



(116) 



where all the digits agree with the numbers obtained for the fit (fewer digits are 
displayed for the results in the table). 

6. 3. Casimir energy of an inhomogeneous string 

The energies (frequencies) of a inhomogeneous string depend both on the 
density of the string itself and on the extension of the string; when the density 
or the length of the string (or both) are changed, these energies change in a way 
which can be described in our WKBPT scheme. In particular we have shown 
that one can extract the asymptotic contributions to the energies of the string 
from the perturbative corrections in the limit n> 1, with a rigorous procedure 
which does not involve any approximation: therefore one can obtain arbitrarily 
accurate expressions for the energies of the string as functions of n working to 
high orders in WKBPT. 

We will now apply this result to the calculation of the Casimir energy of 
an inhomogeneous string. The problem of calculating the Casimir energy of a 
piecewise uniform string has been originally considered by Brevik and Nielsen 
in ref. 11 1 and by later studied in a sequence of works 12, ]]J 14, 15 1. 

Here we consider a string of arbitrary density with Dirichlet boundary con- 
ditions; the Casimir energy of the string is defined in a standard way as (15| 



Ec = 2 E ( Wn ~ w ") 

n=l 



(117) 



where w n are the frequencies the particular string considered and ZJ n are the 
frequencies of the string when its length is made infinity and the density con- 
stant. 
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Working with the energy to order 1 /n 2 we can write explicit expressions for 
the frequencies: 



a{L) 



2mr 



8n 3 7r 3 



a{Lf 



+ 0[l/n 5 l ,(118) 



where the general form of ^3 has been already derived using WKBPT. 

Coming to the calculation of the Casimir energy, we may use a simple cutoff 
function f(n) — e~ an as done in ref. [11] to regularize the divergent behaviour 
of the infinite series: 



C(3) 



- + — + 0[a 2 
12 a z 



<r(Lf 
8tt 3 



47T 2 A 3 



{Vf 



a(L) 



loga + 0[a])^(V) 



2tt 



1 OO 



(119) 



Although it is not wise to attach a physical meaning to this expression until 
a consistent renormalization is carried out (which is beyond the goals of the 
present discussion), leading to a result which is free from divergencies, we may 
observe that the terms stemming from the higher order expansion in powers of 
1/n lead to unambigous finite expressions which are functionals of the density. 
This is the case of the term proportional to C,(3). The divergent behavior (in 
the limit a — > 0) is contained in the first two terms and it should be cancelled 
by the last term. 

6.4-. An asymmetric density 

We now study the example considered by Bender and Orszag at pag. 491 of 
ref. [H, corresponding to the density (we center the string in the origin and set 
L = it/2) 



p(x) 



(120) 



In this case 



8x* + 36nx 2 + 54tt 2 x + 19tt j 
24 



and 



(121) 



V(x)=- 



. r + - r] '' (3cr + 7T 3 ) 2 



3tt> 



(122) 



Lindblom and Robiscoe have also used this example in ref. [6| to test the 
accuracy of their method. 
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Before discussing in detail the application of the results of the previous sec- 
tions to this model, we should obtain accurate approximations to the solutions of 
this inhomogeneous string, which will be needed in order to assess the usefulness 
of our methods. 

A general approach, which can be used for strings of arbitrary density, is 
based on the discretization of the Helmholtz equation on a homogeneous grid 
using a collocation approach based on Little Sine Functions (see for instance 
ref. [l||). Using this approach one is able to obtain very precise values for the 
energies and eigensolutions of a large number of states of the inhomogeneous 
string: in particular, the energy of the fundamental mode of the string with the 
density of eq. (|120[) using a fine grid of 2500 points is 



E, 



LSF 250 




0.001744013543 



(123) 



where all the digits displayed are expected to be exact. Notice that this approach 
is more precise for the lowest states. 

While the collocation approach is general, for the example of a string with 
density (|120[) we may obtain extremely precise results observing that the matrix 



elements of the operator O 



in the basis of the eigenfunc- 



/ p(x) ax p( x ) 

tions of Q can be calculated analytically. This procedure amounts to using a 
spectral method compared with the pseudospectral method used in collocation. 
We report here the results: 



\Ol 



= 5„ 



9n 2 
49tt 4 



4nCi(^)sin(2Hi) 4nCi (^) sin (^) 



49tt 5 



49tt 5 



49tt 5 



49tt 5 



4nSi(2^)cos(2^i) 4nSi cos (2^a) 



for n = m and 



(n\0\m) 



49tt 5 



(to — n) Ci 



. / (to — n)7r 



7 



Ci 



(to — n)n 



, , „(m + n)n\ ({m + n)-K 
+ (m + n) Ci | I — Ci 1 



+ (to — n) Si 



7 

(to — n)7r 
7 



Si 



. , .'(to + n)Ti . 
+ (ro + n)(Sip ^ — ^ — ) S n 



7 



. / (to — n)-7r 



7 

[m + njn 
7 



7 



sm 



cos 



7t(to — n) 



7 



7t(to + n) 
7 

7t(to — n) 
7 

7t(to + n) 
7 



for n =/= to. Here n,m = 1,2,... and Ci(a;), Si(x) are the cosine and sine integrals 
respectively. 

Using these matrix elements we have obtained a N x N matrix representing 
O in the subspace spanned by the lowest N eigenfunctions of Q; the eigenvalues 
and eigenvectors of this matrix are expected to converge to the exact eigenvalues 
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Figure 1: The difference E 



(N) 
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,(400) 



for the lowest eigenvalue of the matrices representing 



OonaJV dimensional subspace. 



and eigenfunctions of O as the dimension of the subspace is increased. This 
behavior is illustrated in Fig. [1] where we have calculated e[ for matrices up 
to size 200 x 200 and we have then plotted the difference e[ N ^ — £^ 400 ^ using 
the very precise numerical eigenvalue of the 400 x 400 matrix. This plot shows 
that the error for the lowest eigenvalue of the 200 x 200 matrix is of order 10~ 18 , 
and therefore much smaller than the one obtained with LSF. 
In this case we have 

E[ 2 00) = 0.00174401354381855 (124) 

where all the digits are expected to be correct. 

This approach allows one to obtain very precise approximations also for the 
excited states of the string: in Fig. [2] we see that the 200 eigenvalues of the 
200 x 200 matrix representing O (circles) are very precise, even for the highest 
states. The few highest eigenvalues in the plot are obtained even more precisely 
using the diagonal matrix elements (n|0|n) (dashed line), corresponding to using 
WKBPT to first order. 

Actually one can use this spectral approach to obtain accurate results also 
for arbitrarily excited states: the crucial observation is that the matrix repre- 
senting O becomes diagonally dominated for highly excited states. This can be 
understood in general terms remembering that the functions <p n (x) of eq. (|23p 
become eigenfunctions of O for n — > oo . 
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Figure 2: The difference E n — E„ for the lowest 200 eigenvalues of the string of density 
(11201) . The circles correspond to E n calculated taking the eigenvalues of the 200 X 200 matrix 
representing O, whereas the triangles correspond to E n obtained using LSF with a grid of 
2500 points. The dashed curve corresponds to the difference |(ti|0|?i} — E n ^|. 



Therefore, if we build a square (2n + 1) x (2n + 1) matrix, centered on 
the n state, its eigenvalues and eigenvectors will provide increasingly better 
approximations to the exact eigenvalues and eigenfunctions of O as n goes to 
infinity. In this way one may obtain precise estimates of the asymptotic behavior 
of the energies of this problem for n — > oo. 

For example, we may consider the three states corresponding to n = 10° — 1, 
10 5 and 10 5 + 1 and build a matrix centered at each of these states with n = 100. 

To calculate the eigenvalues we have worked with the numerical matrices 
with a precision of 100 digits and we have then extracted the coefficients of the 
asymptotic expansion with a fit, obtaining 

E n k 0.0018855885773910003409 n 2 - 0.00026004036832396307366 

0.00075459566769573983984 
+ ^ + ... . (125 

These numerical coefficients should be compared with the exact coefficients 
calculated with WKBPT: 

E n = Am 2 + A 2 + ^ + ... (126) 

where 

Ax = —-T « 0.0018855885773910003409 (127) 

49 7T 
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A, = - 



4tt 6 



-0.00026004036832396307402 



and 



4 1} +4 2) . 



(128) 



(129) 



The contribution to A 3 stemming from first order is easily calculated and 
reads: 



4(1) _< L ) 
A3 ~ ^ 



3577 



V {1) {a{L)) -V (1) (0) = 0.000736 . (130) 

J 5127T 8 



As we have seen before the contribution to A 3 stemming from second order 
involves a series 



(2) a(L) 2 1 



T 4 ^ 2k 2 

k=l 



- T 



2irk 



,a{L))T 



<y(L) 

irk 
2nk ' 



ML) 



<L) 
irk 



5>fc > 



(131) 



k=l 



and requires the calculation of T . 
For this example we find 



G(v,a ) = 



12 (3(To + 



(9 (tTg + u 2 ) + 67r 3 CT + 7r 6 ) 



(132) 



and 



-2ACi(A) sin(A) + 2ASi(A) cos(A) -ttAcos(A) + 2 , s 
F(v,<to) = — o-r (133) 



9u 2 A 



where 



A = 



7T 3 + 3CT 

3v 



(134) 



After substituting the expression for T inside the expression for A^ we 
obtain an explicit expression for jk- 



Ik 



141127T 8 

87rA;Si 



( 

fc-7r\ 
~7~ 



87rfcCi 



cos 



kn 

Y 

nk 

T 



sin 



nk 

T 



+ 87rfcCi 



- 87rfcSi 



. f8kw 
1 — 



cos 

2 



8fc7T 



sin 



a 2, f^ k \ A 2, /87Tfc\ „, 

+ 47T fc COS ( — J + 47T fc cos ( J - 6c 

- 32 ( 2nkCi i— sin (— - 2nkSi i— cos f — 
+ 7r 2 fccos I — J - 7 J • I 167r/cCi I j sin I — ^— j 



- 167rfcSi 



16fc7T 



COS 



167rfc 



8-7r 2 fc cos 



167rfc 



- 7 
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(2) 

Wc may estimate A3 by restricting the series in ' to the hrst N terms. 
For example, using TV = 8000 we obtain 

„___ 8000 

A 3 = A ( } } + 4 2) w — — - + V 7fc w 0.00075459642403370616977 . (135) 

5127T 8 z — ' 

fe=l 

A much better estimate may be obtained taking into account the asymptotic 
behavior of the 7^. for k — > 00: 



and therefore wc may write 



631561441 

7fc ~ 294912^1^4 ( 136 ) 



A (2) _ 



(137) 



^ _ 631561441 ^ r 631561441 

' 2-^ lk ~ 26542080tt 8 + ^ 7fc ~ 2949127r 12 fc 4 
fc=i fe=i L 

If we limit the sum to the first 8000 terms as before we obtain a more precise 
estimate for A 3 : 

(i) 631561441 y!° I" 631561441 " 

3 ~ 3 + 26542080tt 8 + ^ lk ~ 2949127r 12 fc 4 

fe=l L J 

w 0.00075459642403521434583 . (138) 



Therefore we see that the coefficients obtained with the fit are in very good 
agreement with the theoretical values obtained with WKBPT. 

We will now use the very precise numerical results at our disposal to dis- 
cuss the WKB-perturbation method and its improved version, which we have 
described in the previous sections. 

We first concentrate on the bounds for the energy of the fundamental mode, 
starting with the one obtained using DPT, which largely overestimates the exact 
result 

Ex < — Si(27r) - — Si(47r) + — — r + — « 0.002868007277 . (139) 

157T 15-7T 307T 4 157T^ 

We should not be surprised by this result since the density of the string 
varies strongly on the domain and therefore density perturbation theory is not 
applicable. 

We now discuss the first bound obtained WKBPT. First we make an im- 
portant observation: since V(a) is negative on the whole string, then E„ in 
eq. (|59j) is also negative. Moreover, for the modes of the string which are arbi- 
trarily excited the ratio / becomes arbitrarily small, while E^ becomes 
arbitrarily accurate as n — > 00, thus we may conclude that the E^ are upper 
bounds to the energies of the string as n — > 00. We have numerically tested this 
conjecture using the eigenvalues of the 400 x 400 matrix representing O and we 
have found that it actually holds even for the intermediate states. 
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The WKBPT bound is readily obtained using the matrix element previously 
calculated: 



El < »* 



49tt 5 49tt 5 



Cil^-Cif^W^ 



7 \7 V 14 I 



2tT\ „. /167T\\ . /37T 



Si — - S 



7 \ 7 V 14 



sin 



0.001745807766 , (140) 



which is much tighter than the bound obtained with DPT. It is also interesting 
to compare this result with the WKB result to first order: 

£ WKB(0) ^WKB(l) = 9 _ 1 w 0Q162555 (W1) 

1 1 49?r 4 4tt 6 v ; 

which is less precise than the previous one and falls below the exact value. 

As we have discussed before, one may still obtain better bounds using the 
first order WKBPT wave function as trial solution: in the present example, the 
matrix elements (A;|V|i) are known explicitly and therefore the bound can be 
calculated easily. Notice that the variational bound holds also when a finite 
number of terms is used in the series appearing in (O)*. 

In Table [5] we report the variational bounds obtained restricting the sum up 
to k = N terms and compare the results with the "exact" numerical value quoted 
by Bender and Orszag and with the even more precise values that have been 
obtained using a collocation approach based on the LSF functions of ref. [l6| 
with grids with 2000 and 2500 points (the use of collocation in the numerical 
solution of the problem of an inhomogeneous string is illustrated in ref. [l[). 
Notice that for N > 20 the variational bound does not improve, signaling the 
need for a more refined ansatz. 

In Table |3] we report the energies of selected modes of the string considered 
by Bender and Orszag, calculated with WKBPT to first order and second orders 
(second and fourth columns), with the asymtotic formulas Ain 2 +A2 and Ain 2 + 
A-2 + (third and fifth columns) and with the LSF collocation method with a 
grid with 2500 points (sixth column) . Clearly the second order WKBPT results 
are very precise even for the lowest modes: this fact signals that the dominant 
contributions to the asymptotic coefficients Aj, with j — 4, 5, . . ., may be (at 
least in principle) calculated using the expression of the energy to second order 
in WKBPT, in analogy with what we have done in the calculation of A3. 

We will now discuss the bound obtained using the improved WKBPT ap- 
proach (iWKBPT), eq. (|7T|) . We consider the trial density 



p{x) = p{x) 



N 



(142) 



where the a,j are N + 1 real parameters to be determined variationally. 

Clearly choosing a = 1 and dj = for j > (corresponding to N = 0), 
p{x) reduces to the physical density: in this case we recover the simpler bound 
of WKBPT. 
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Table 2: Variational bounds using first order WKBPT to the energy of the fundamental mode 
of a string with density p(x) = (x + ^7r) . The notation [— n] means 10 _ " . 



N 


Et(N) 


Ex(N) 






2 


0.0017442945174961415 


2 


81 


[-7] 


3 


0.0017440904942458490 


7 


70 


[-8] 


4 


0.0017440396134339554 


2 


61 


[-8] 


5 


0.0017440245258144354 


1 


10 


[-8] 


6 


0.0017440189344824756 


5 


39 


[-9] 


7 


0.0017440167007034352 


3 


16 


[-9] 


8 


0.0017440156851693774 


2 


14 


[-9] 


9 


0.0017440152042467934 


1 


66 


[-9] 


10 


0.0017440149552016950 


1 


41 


[-9] 


20 


0.0017440146360548171 


1 


09 


[-9] 


30 


0.0017440146309786759 


1 


09 


[-9] 


40 


0.0017440146306353234 


1 


09 


[-9] 


50 


0.0017440146305891107 


1 


09 


[-9] 


ref. [3] 


0.00174401 




LSF2000 


0.0017440135432079554 








LSF2500 


0.0017440135430055502 
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Tabic 3: Comparison of the energies for selected modes of a string with density p(x) = 
(x + | tt) calculated with WKBPT to first and second orders (second and fourth columns) , 
with the asymtotic formulas Ain 2 + A2 and A\n 2 + A2 + ^ (third and fifth columns) and 
with the LSF collocation method with a grid with 2500 points (sixth column). The series in 
the second order WKBPT term is restricted to the 20 closest states. 



n 


J 


Am 2 + A 2 


IT 11 


Ain 


2 _l_ A J- A 3 


77'LSF2500 


1 





00174581 





00162555 


0.00174405 





00238014 





00174401 


2 





00734927 





00728231 


0.00734864 





00747096 





00734866 


3 





01675247 





01671026 


0.01675237 





01679410 





01675238 


4 





02993818 





02990938 


0.02993827 





02995654 





02993828 


5 





04690044 





04687967 


0.04690060 





04690986 





04690060 


6 





06763676 





06762115 


0.06763693 





06764211 





06763693 


7 





09214593 





09213380 


0.09214609 





09214920 





09214609 


8 





12042730 





12041763 


0.12042744 





12042942 





12042744 


9 





15248051 





15247263 


0.15248064 





15248195 





15248064 


10 





18830534 





18829882 


0.18830546 





18830636 





18830546 


20 





75397717 





75397539 


0.75397721 





75397728 





75397721 


30 


1 


69677048 


1 


69676968 


1.69677050 


1 


69677052 


1 


69677050 


40 


3 


01668214 


3 


01668168 


3.01668215 


3 


01668216 


3 


01668215 


50 


4 


71371167 


4 


71371140 


4.71371170 


4 


71371170 


4 


71371170 
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Table 4: The variational estimate for the energy of the fundamental mode of the string of 
density p(x) = (x §•?>") 4 and the quantity ^p(x) — ^ x determining the asymptotic 
behavior of the spectrum. 



iV 


^iWKBPT 


^iWKBPT 


_ ^(200) 


a(L) 2 r 
aiL)* J- 


-L p(x) aX 





0.00174580776578874 


1.79 


[-6] 


1 


00000 


1 


0.00174408029127495 


6.67 


[-8] 


1 


00076 


2 


0.00174402710509883 


1.36 


[-8] 


1 


00088 


3 


0.00174401553013316 


1.99 


[-9] 


1 


00088 


4 


0.00174401377197412 


2.28 


-10] 


1 


00088 


5 


0.00174401354435309 


5.35 


-13] 


1 


00088 


6 


0.00174401354405299 


2.34 


-13] 


1 


00088 


7 


0.00174401354384465 


2.61 


-14] 


1 


00088 


8 


0.00174401354382102 


2.47 


-15] 


1 


00088 


9 


0.00174401354381886 


3.12 


-16] 


1 


00088 


10 


0.00174401354381859 


3.70 


-17] 


1 


00088 


£,(200) 


0.00174401354381855 







In Table |4] we report the energy corresponding to a set of values of the 
"optimal" aj and the numerical value of the quantity ^Tgp S-l ^p(x) wn i cn 
determines the asymptotic behavior of the energies. 

For N = 10 the first 17 digits of the variational energy agree with those of 
the energy of eq. (|124|) ; the case N — corresponds to normal WKBPT, which 
uses the physical density as p{x). 

In Fig. [3] we plot the physical density p(x) = (a; + |7r) (solid line) and 
compare it with the optimal trial density obtained with the iWKBPT bound 
using N = 10 parameters (dashed line). This density approximates the effective 
density for the fundamental mode of O. 

Using the transformation discussed by in ref. 0] we may obtain refined 
bounds for the energy of the fundamental mode of our string, studying an 
isospectral string of density 

16(tto + l) 2 (^(2a(2x + tt) + 3) + 2x) 4 

P\ x ) = i—R — ; — tt^s ■ ( 143 ) 

[a(2x + it) + 2) s 

In Table [S] we report the values of these bounds obtained by minimizing the 
expectation value of O in the fundamental mode of Q where the density p(x) 
has the same form considered before. The minimization is done with respect 
to a and to the parameters in p(x). In the fifth column we report the optimal 
value of a, for a set of N parameters a^. The minimization of the expectation 



34 



1200 



1000 




1.0 



Figure 3: Comparison between the physical density p(x) = (x + |tt) (solid line) and the 
optimal trial density obtained with the iWKBPT bound with N = 10 (dashed line). 



value for a set of N + 1 parameters (the N parameters a, and a) is done using 
the N parameters of the previous solution as a starting point. 

Notice that the result obtained for N = is roughly an order of magnitude 
more precise than the corresponding result in Table |4] and that it retains the 
exact asymptotic behavior. 

On the other hand the results obtained for N > 7 have a larger error than 
the corresponding results obtained in Table [4j this means that the solution that 
we have found corresponds to a local minimum of the expectation value. On the 
other hand, we should observe that the asymptotic behavior of the spectrum 
obtained in table [5] for N > 7 more accurate than the corresponding behavior 
in Table |U 



6.5. A string with rapidly oscillating density 

We will now study a string of unit length (|x| < 1/2) and with a rapidly 
oscillating density 

>2 



p{x) = (2 + sin(1007ra))' 



In this case one has 



a(x) = 1 + 2x + 



sin (50tt:e) 
50tt 



(144) 



(145) 



It is known that the solutions for strings with rapidly oscillating density have 
peculiar properties: in particular the solutions whose wavelength corresponds 
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Figure 4: \E l ^ VKBPT - E^ xact \ for the case of a string with density p(x) = (x+ §7 
iWKBPT results are obtained using an optimal density with N = 10 parameters. 
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Table 5: The variational estimate for the energy of the fundamental mode of the string of den- 
sity p(x) = m ^ a+1 \^f"^l)^2)« +i)+2X) and the q uantit y f^jr S-L P{ %) dx determining 
the asymptotic behaviour of the spectrum. 



N 


^iWKBPT 


^iWKBPT 


_ ^(200) 


*(Lf r+L p(xfl^ 
a(L) 3 J-L p(x) 


a 





0.00174418771063455 


1.74 


[-7] 


1.00000 


-0.00840 


1 


0.00174404317510904 


2.96 


[-8] 


1.00415 


0.01082 


2 


0.00174401638359089 


2.84 


[-9] 


1.00532 


-0.02892 


3 


0.00174401358072024 


3.69 


-11] 


1.00930 


-0.03528 


4 


0.00174401356813542 


2.43 


-11] 


1.02434 


-0.05076 


5 


0.00174401354405135 


2.33 


-13] 


1.00016 


-0.00565 


6 


0.00174401354405132 


2.33 


-13] 


1.00016 


-0.00565 


7 


0.00174401354385037 


3.18 


-14] 


1.00016 


-0.00566 


8 


0.00174401354382192 


3.37 


-15] 


1.00016 


-0.00566 


9 


0.00174401354381896 


4.10 


-16] 


1.00016 


-0.00566 


10 


0.00174401354381860 


5.05 


-17] 


1.00016 


-0.00566 


^(200) 


0.00174401354381855 
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Figure 5: Wave numbers of the string with density p(x) = (2 + sin(1007r:r)) 2 obtained with 
LSF with a grid of 2500 points. The dashed line is the asymptotic behavior k n = nir/a(L). 



roughly to the typ ical size of the density oscillations are localized at the ends of 
the strings [3. Il7l Hc|. 

Before we apply the results of the previous sections to the study of this 
string we wish to obtain numerical results which will be useful to assess the 
precision of the analytical methods: as for the previous example we have used 
a collocation approach based on LSF with an homogeneous grid of 2500 points. 
The grid spacing h = 2L/N is sufficiently fine to obtain precise results for the 
first few hundreds of states; in particular, the energy of the fundamental mode 
obtained with this grid is 

e lsf 25QQ = 2.1931584082 . (146) 

The quality of the variational bounds obtained with the different approaches 
will be established using this value. 

In Fig. [5] we display the numerical values of the square roots of the first 350 
energies of this string obtained with the LSF collocation using 2500 points (plus 
symbols) and we compare them with the exact leading asymptotic behavior 
k n = mr/2 (dashed line). Gaps are observed in correspondence to quantum 
numbers which are multiples of 50: as we have anticipated these gaps occurr 
because of the appearance of localized solutions at one of the ends of the string. 

In Fig. [S] we display three solutions of the Helmholtz equation for the string 
with density p{x) = (2 + sin(1007r:r)) 2 , obtained using LSF collocation with a 
grid of 2500 points. The mode for n — 50 is localized on the left end of the 
string, while the modes corresponding to n — 49 and n — 51 extend over all the 
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string. It is instructive to decompose these solutions in terms of the modes of 
an uniform string: we may expect that the modes of this string which are not 
localized may be described in terms of few modes of an uniform string, whereas 
the modes of this string which are localized necessarily involve a great number 
of modes of the uniform string. 

This analysis can be carried out rapidly and efficiently using the results of 
the LSF collocation, and at no extra computational cost. The fundamental 
observation is that the LSF are defined as 

2L N 

Sk(x) = — 22 tpn{xk)lpn{x) , (147) 
n=l 

where i[> n (x) = y^"sin n7! ^+ L } are the eigensolutions of the uniform string on 
— L < x < L with Dirichlet boundary conditions at x = ±L. The Xk = 2Lk/N 
are the uniformly distributed grid points. 

A function fix) whose values are known at the grid points may then be ap- 
proximated at arbitrary points in the domain through the interpolation formula 

N/2-1 

f(x) ss ^2 f{xk)sk{x) 

fc=-7V/2+l 

AT/2-1 



N 

E 

n=l 



Jj E f( x h)^n{Xk) 
fc=-JV/2+l 



M*) > ( 148 ) 



where the terms inside the parenthesis are the approximate Fourier coefficients 
of f(x). In the case that we are considering f(xk) is just the k th entry of one of 
the eigenvectors of the (N — 1) x (N — 1) matrix that we have obtained from the 
discretization of the Helmholtz equation, multiplied by a factor yj N/2L, and 
ipn( x k) is the k th entry of the vector obtained evaluating the n th eigensolution 
of the uniform string on the k th grid point. The calculation of the Fourier 
coefficients thus involves (apart from a multiplicative factor) the multiplication 
of two vectors, which have been already calculated. 

In Fig. [7] and [8] we display the approximate Fourier coefficients for the so- 
lutions of Fig|6] for n odd and even respectively, obtained with LSF on a grid 
of 2500 points. The coefficients of the solutions corresponding to n = 49 and 
n = 51 behave very differently for odd and even modes while the coefficients 
of the localized solution, corresponding to n — 50, behave similarly in the two 
cases and decay slowly. 

We may now discuss the bounds for the energy of the fundamental mode 
obtained using the different approaches outlined in the previous sections; in the 
case of DPT one uses the fundamental solution of a uniform string as variational 
ansatz and obtains a bound which falls quite above the exact energy: 

E x < 6335.15 . (149) 
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Figure 6: Solutions for the modes 49,50 and 51 of the string with density p(x) 
(2 + sin(1007nz)) 2 obtained with LSF with a grid of 2500 points. 
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Figure 7: Approximate odd Fourier coefficients for the modes 49,50 and 51 of the string with 
density p(x) = (2 + sin(1007rz)) 2 obtained with LSF with a grid of 2500 points. 
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Figure 8: Approximate even Fourier coefficients for the modes 49,50 and 51 of the string with 
density p(x) = (2 + sin(1007rx)) 2 obtained with LSF with a grid of 2500 points. 



The situation is only partially better using WKBPT: the "potential" V{x) 
corresponding to this density is 



V(x) = - 



12507T 2 (8sin(1007rx) + cos(200?ra;) + 5) 



(sin(1007re) + 2) 4 
and using the formula derived earlier the bound reads 

E 1 < 1189.60 . 



(150) 



(151) 



This large bound originates almost entirely from the first order contribution 

L / \ i t-L 



JpJx)V{x)W^dx 
v{L) J_ L a{L) 



1 



y r p(x)V(x)dx . (152) 



This result clearly suggests that the WKBPT describes poorly the lowest 
part of the spectrum of this system, although the basis of the eigenfunctions of 
Q reproduces the asymptotic behavior of the high part of the spectrum of O. 
Both bounds calculated so far are clearly useless! 

In the iWKBPT approach one may use an effective density p(x), different 
from the physical density and specific to a given state, to improve the descrip- 
tion of the lowest energy solutions. As we have seen in the case of WKBPT 
discussed above, the large values obtained in the bound originate from the rapid 
oscillations of the density: one should therefore pick an effective density p(x) 
for which these contributions arc suppressed. 
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An useful ansatz is 



p(x) w p(xf 



(153) 



since it cancels large contributions to the expectation value of O stemming from 
the oscillations of the density. 

The bound for the fundamental energy obtained in this case is much closer 
to the numerical result obtained with collocation 



#i < 3.07325 



(154) 



although the quality of this bound is not competing with the precision of the 
numerical result. 

We will now show that it is possible to obtain extremely precise variational 
bounds, choosing a suitable variational ansatz, given by the trial function (nor- 
malized to one) as 



V2 



\J p(x)/L si 



sm 



L 



2L 



(155) 



A simple calculation shows that the bound is this case is 



E X < — 7T 

9 



2.19325 , 



(156) 



where the first four digits of the rhs are correct. 

Can we still improve this bound? Yes, we can: we can use Theorem 1 of 
ref. [l[ to obtain further analytical improvements of this bound. Using the nota- 
tion of that paper we set £o(%) = ipi {x) and obtain refined approximations 
for the fundamental solution applying the iterative equation eq. (|92|) . 

The bounds obtained in this way are given in Table |H1 after 7 iterations it 
appears that the first 19 digits have converged. In principle the iterations can 
be performed to arbitrary orders, as long as one is able to perform the integrals 
analytically 0. A comparison with the LSF result shows that the first 10 digits 
of the value obtained with collocation are correct. 

Having obtained accurate numerical solutions for the low lying states of this 
string, we may look for the effective densities corresponding to these states. In 
this case one needs to solve numerically the equation 



(j{x) 1 . 2t:5{x) 

J{L) = 2> Sm a{L) 



-L 



i (trial) / \ 



dy 



(157) 



The effective density is then obtained through the relation a'(x) — J pix). In 
Fig(3] we show the effective densities for the first four solutions, scaled by the 
factor p(x) 2 , which we have already seen to be a good approximation: as a 



7 The form of the density in this example was chosen to allow to perform several iterations. 
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Table 6: Variational bounds for the energy of the fundamental mode of the string of unit 
length of density p{x) = (2 + sin(1007rx)) 2 obtained using Theorem 1 of rcf. [ij with the 
initial trial function in eq. 1 11551 . 



n 


E x 





2.1932454224643019153 


1 


2.1931584087453807530 


2 


2.1931584087440065685 


3 


2.1931584087439513989 


4 


2.1931584087439479612 


5 


2.1931584087439477463 


6 


2.1931584087439477329 


7 


2.1931584087439477320 



matter of fact we may notice that p\(x)j p 2 (x) is almost fiat, with oscillations 
only at the ends of the string. 

In the case of the remaining states we notice a very peculiar behavior of the 
effective density, which develops n—1 zeroes in the domain, which also coincide 

with the zeroes of the solution (the zeroes of the eigenfunctions of Q correspond 
either to the zeroes of the trigonometric function or to the zeroes of the effective 
density, as in the present case). 

In Fig. [10] we display the same ratio for the modes corresponding to n = 49, 
50 and 51. The effective density of the mode 50 is clearly localized at the left 
end of the string. 



7. Low and high energy asymptotics for strings with rapidly oscillat- 
ing density 

In this section we use the theorems of Section [5] to obtain the low energy 
asymptotics of a string with density § 

p(x) = 2 + sin V - 2 > (158) 

with e -> 0+ and x 2 < 1/4. 

As we have seen through Theorem [JJ it is possible to obtain a sequence of 
functions which converge to the fundamental solution of O provided that the 



8 The low energy asymptotics for this problem has been obtained by Castro and Zuazua in 

a. 
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Figure 10: Effective densities for the solutions 49, 50 and 51 of the string with density p(x) = 
(2 + sin(1007nz)) 2 obtained with LSF with a grid of 2500 points. 
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initial trial function is not orthogonal to it. Therefore we may consider the 
function 



= v/pRshr n{X ^ /2) • (159) 

Actually this £o(x) not only is not orthogonal to the exact fundamental 
mode of O, but it also tends to become an exact eigenfunction of O in the limit 
0+. 

We can see this by calculating 



having substituted in the last expression the highly oscillatory term sin 



(160) 



with its average value on the interval. 

We thus use Theorem [T] to build increasingly refined approximations of the 
fundamental solution of O corresponding to the density (I158|) in the limit e — > 
+ . For instance, after one iteration we obtain the function: 



2tt 2 (e 2 - 4) 2 



4e 3 (2x + l)cos 



+ e 2 ^4e(2:r-l) + (e-2) 2 sin 
+ (e + 2) 2 sinl -))+4(e -4) co$,{itx) (161) 

We have also calculated explicitly the functions corresponding to two and 
three iterations, although their expressions are much lengthier and therefore we 
will not report them here. 

We will instead calculate the expectation value of O in these functions, thus 
obtaining an asymptotic expression for the fundamental energy of this string in 
the limit e ->• 0+. 

We report here these expressions: 

,Ax f-L €a{x)0€o(x)dx 7T 2 1 . 2 /7T\ 3 1 . 2 /7T\ 5 . 

<°>* s tf^y* " T + 4— (l) ^ + (l) g5 

V /6 ^ii{x?dx 2 64 4 U 1024 

7r(5sin(^L) - 116 cos (^) + 116) - , s 

+ Kt > 1Q24 Kt > L ^ + --- (163) 
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(0)i 



S-L Mx)dW)dx _ IT 2 7T 2 2 

£Z&(x)*dx * 2 64 6 
7T (5 sin (^) - 116 cos (^) + 116 
1024 



— 7r sin — e 



15tt 2 
1024 ' 



(164) 
(165) 



As we can see the expectation value of O has converged to order e 5 after 
just one iteration; thus we can conclude that, for e — > + , the energy of the 
fundamental mode goes like 



7T 2 7T 2 9 1 , /7T\ 

Ei - T-6i e + i™ (7) 



e 3 - 



15tt 2 
1024 ' 



(5sin(^F) - 116 cos ^f) + 116 
1024 



0[e«] 



(166) 



This result may be compared with eq.(3.33) of Ref.0], which reports \fE\ 
to order O [e 4 ] for this example. 

Let us now derive the asymptotics for the low energy states of this string; 
in this case we use Theorem [3] and consider the trial function 



Vo(x) = 



\ 



2 " (« + j: 



7r a; 



sin 



+ !)' 



with rjo(-L) = f]o(L) = 0. Notice that r]o(L) = only for L 
n = 1, 2, . . . and 

6770(2;) n 2 ir 2 n 2 ir 2 



mix) 



(167) 



— , with 



(168) 



meaning that for e — > + , the 770(2;) for L 



- converges to the n th solution 



for this string (again we are substituting the oscillating factor with its average 
value on the interval in this limit). 

We may now apply Theorem [3J after one iteration we obtain the function 



mix) 



(2L + l) 2 Jsin 



27r(a;+i) 



87r 2 (e-2L-l) 2 (e + 2L + l) 2 



4e 3 (2a; + 1) cos 



/2ttL + tt 



e 2 (e + 2L + l) 2 cos 
e 2 (e- 2L- l) 2 cos 



n(2x + l)(e - 2l - 1) 
2eL + e 

vr(Z(-2e + 4x + 2) + 2(e + l)x + 1) 



2eL 



+ 4 ( 2e 3 (x — L) + (e 2 — (2L + l) 2 ) 2 sin 



2irx + 7r 
2L + 1 



(169) 
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Figure 11: The function of eq. l|169|l evaluated at x = L, for e = 1/50 as a function of (L + L). 



The zeroes of the equation 



m(L) = o 



(170) 



will now provide improved values for L and for the energies. 

In FigfTT]we display rji(L) corresponding to e = 1/50, as a function of (L+L): 
the spikes of the curve correspond to zeroes of rji (L), located at specific values of 
L. Corresponding to each value Theorem 3 provides an approximate eigenvalue 
and eigensolution to the string problem. 

In Table [7] we report the values of L corresponding to the first 20 zeroes of 
7/1 (L) and the corresponding energies obtained by calculating the expectation 
value of O in the corresponding state; the last column contains the precise 
numerical results obtained with LSF collocation with a grid of 2500 points. 
Thus we see that Theorem 3 provides rather precise results for a large number 
of states just after one iteration. 

We look for approximate solutions of the form 



L 



1 1 

2 + n 



SL 



(171) 



where SL is a correction. To leading order we have 

e 3 (-l)-"n (racos(2|) + (-1)" cos (2s) 



SL 



8tt 



n - 1 



(172) 



Notice that SL = for n — 1, as it should be. 
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Table 7: Approximate energies of a string with density H158H with t = 1/50 obtained with 
Theorem [3] (third column) and corresponding values obtained with LSF on a grid of 2500 
points (fourth column). The second column contains the values of L for which r]i(L) = (see 
the figure). 



n 


L 


E n 




1 


0.50000000000000 


4.93474049237786 


4.93474049164385 


2 


0.00000254847810 


19.7382203568987 


19.7382203589595 


3 


-0 


16666618871377 


44.4082063869071 


44.4082063970937 


4 


-0 


25000000000000 


78.9409479076410 


78.9409485235574 


5 


-0 


30001279890522 


123.331046653456 


123.331137037346 


6 


-0 


33332758093170 


177.571827239743 


177.571841516164 


7 


-0 


35716402051048 


241.652567358357 


241.654430048229 


8 


-0 


37497950785717 


315.566145317546 


315.568467036379 


9 


-0 


38888465796479 


399.301373439805 


399.301587271749 


10 


-0 


39993517619554 


492.724112367784 


492.839343430514 


11 


-0 


40908535234753 


596.163807132223 


596.165023435926 


12 


-0 


41659695736991 


708.868637232727 


709.259433238896 


13 


-0 


42315478442053 


831.113631633346 


832.100639666984 


14 


-0 


42856494498590 


964.660849364111 


964.663666640351 


15 


-0 


43336143492296 


1106.57429450867 


1106.92013667688 


16 


-0 


43735701119023 


1249.09849824817 


1258.83784766743 


17 


-0 


44135393715531 


1396.14689996551 


1420.38027258297 


18 


-0 


44431129182849 


1574.93187010870 


1591.50596694241 


19 


-0 


44739462363259 


1770.83890618558 


1772.16786511172 


20 


-0 


45000000000000 


1962.04796373469 


1962.31244196980 
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Working to order e 3 we have calculated the expectation value of O in 771 (x) 
to order e 3 : 

<0> " = "l^Si^ - ~ _a " 7 (7) H 

which provides the asymptotic behavior of the energy of this string for e — > + . 
Notice that 

/— 7m 7re 2 n 3 e 3 n 3 sin 2 (f) 

coincides with the result obtained by Castro and Zuazua in ref. [2J using the 
WKB method. The advantage of the present approach is that the calculation 
of the higher order asymptotics only requires to apply further iterations and 
expand the expectation value of O to higher order (of course one needs to be 
able to perform the iterations explicitly!). Moreover, the theorem [3] allows to 
obtain at the same time approximations for the energies and for the solutions. 

We will now discuss the high energy asymptotics for this string, using the 
asymptotic expansion obtained with WKBPT. A simple calculation shows that 
in the present case: 



^ (l 6 sin (22EEfcr) + CQS + 9) 



V(x) = * i 5 — '- . (173) 

8e 2 (sin (*™pL\ +2 ) 3 



We also find that 



V3e (E (f 



2\ _ zp ( 7r(e-4x-2) i 2 
3/ \ 4e 1 3 



a(x) = / y/fajdy = ^ ^ >-L , (174) 



— L " 

where E(x\m) = \J 1 — m sin 2 t dt is the incomplete elliptic integral of second 
kind. 

Notice that 

lim a(L) = — E ( ^] « 1.39066 . (175) 

e->0+ 7T \S J 

The dominant contribution to the high energy asymptotics of this string is 
readily obtained using these expressions: 

7TTL 

E n « ^TT + • • • ~ 2 - 26 « + • • • > ( 176 ) 

for n — > oo. 

We come now to the subleading term, which requires the calculation of 



1 



r(L) 
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23. | 7rdn(g) h ^Qg^j) 

72V2 e( x(L) + l8ea(L) (sin (2s) + 2) 3/2 + 6V3ea(L) 

^(gjgj) ^(fll) 23,cos(f) 

3V3ea(L) 3V3ea(L) 36efT ( L ) ( sin (2l) + 2 ) 3/2 

ttF (esc" 1 (V3) 1 1) 



6V2ecr(L) 



(177) 



where -F(a;|m) = £f , 1 (it is the incomplete elliptic integral of first kind. 

' W l—m sin 2 t 

Notice that e 2 (V) tends to a constant as £ -4 + : lim c ^ o +(^0 = K - There- 
fore, for a given e, arbitrarily small and positive, the behavior of the energy will 
be dominated by the leading term only for states with n 3> ncr(L)/Tre 2 . 

*»e»<V> = £ ( 2 -|||^ 0.214515, (178) 

where K(m) and E(m) are the complete elliptic integrals of first and second 
kind respectively. 



8. Conclusions 

In this paper we have devised an alternative perturbative approach to the 
problem of finding the normal modes of a string with arbitrary density: using 
this approach we have obtained an asymptotic high energy expansion, and we 
have derived explicit analytical expression for the first three coefficients of this 
expansion. The first two coefficients were already known and they can be derived 
using the WKB method; the third coefficient was not known (to the best of our 
knowledge) and it is expressed in terms of a series whose summands are themself 
series (divergent ones), which can be evaluated using a Borel transform. Higher 
order coefficients could also be calculated in a similar fashion, although we have 
not done it here. We have discussed specific cases, where exact solutions (or 
very precise numerical solutions) were available, and we have reproduced the 
theoretical values obtained with our formula with high accuracy. 

We have also applied the iterative theorems of ref. [l| to the problem of a 
string with highly oscillatory density; in this way we have obtained an explicit 
analytical formula which describes the low energy asymptotics of a string with 
density of arbitrarily small wavelength and we have thus reproduced the results 
of Castro and Zuazua 0, who had used the WKB method. 

We believe that several issues could be investigated using the techniques 
described in this paper and in ref.pj; in particular, we plan to calculate in a 
future work few more higher order asymptotic coefficients for the energies of 
strings with arbitrary density using the WKBPT method of the present paper 
and to perform a similar calculation of the asymptotic form of the solutions, 
which has not been carried out here. 
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Finally, we would like to bring the attention of the reader to a " pedagogical" 
aspect of the perturbative scheme devised in the present article: while both 
perturbation theory and the WKB method are a standard part of most quantum 
mechanics/mathematical physics books, normally they are discussed separately 
as unrelated topics. The WKB-pcrturbation method that we have described 
here shows that it is possible to merge the two approaches in a single and more 
powerful method, at least for the problem of calculating the normal modes of 
inhomogeneous strings of arbitrary density. Its extension to higher dimensional 
problems is certainly not trivial, although it deserves to be investigated in the 
future. 



Appendix A. Asymptotic expressions for the the matrix elements 



In this appendix we work out the asymptotic expressions for the matrix 
elements (n|y|fc), which appear in the WKBPT perturbative series. 

We start with the diagonal terms, which appear in the first order correction 



(n\V\n) 



a(L) 
(V)- 



a(L) 



sin 



a{L) 



"TFT / cos 

v( L ) Jo 



V{x(a)) da 
( 2nira 



V(x{a)) da (A.l) 



where 



(V) 



r(L) 



V(x(a)) da 



is the average of V on the interval (0, a(L)). 

The second term may be evaluated using the integration by part: 



(II) 



1 



a(L) 



<x(L) 



COS 



2mra 



V(x(a)) da 



2mr 

A-K 2 n 2 
a(L) 



1 r (L) . /2nira\ dV(x(a)) 
— J Sm {-m)^— (h 



dV{x{a)) 



da 

<x(L) 

COS 



da 
dV(x(a)) 



da 



o 



a{L) 

/2nira\ d 2 V(x{a)) 
(ML)) da 2 



da 



(A.2) 



(A.3) 



Since the last term in this expression has a form which is similar to that of 
the expression that we started with, the integration may be obtained with no 
effort: 



(II) 



= "TFT / cos 
Jo 



(2mra 



V{x{a)) da 



■50 



E^S^" 1 )* V^(*(L))-V^\0) . (A.4) 



fe=0 



(27m) 



Therefore 



^H^-ESm(-i) fc ^ +1) Ki))-^ +1) (0) ■ (A.5) 



fc=0 



(2-Kn) 



The same strategy may be followed for the non-diagonal matrix elements; 
after performing repeated integrations by part one is left with the expression: 



(n\V\k) = 



a(L) 



a(L) 

1 



nna 



a{L) 



a{L) 
(k — n)ira 

<y(L) 



kna 



V(x(a)) da 
(k + n)na 



; a{L) 



2j+l 



> \^(fc_ n )2j+2 ( fc + n )2j+2 

. [(_l)fe+n F (2i+l) ((t(jL)) _ V (2j + 1) (Q) 



V(x(a))da 



(A.6) 



Notice that: 



• lim r , 



lim,, 



,(n\V\n) = (V); 
>{n\V\k) = for k ^ n; 



As a result, E„ dominates all the remaining terms in the perturbative series 
as n — > oo, as it should be. In particular, for u» 1 the non-diagonal matrix 
element behaves as 



(n\V\k) 



cr(L) Akn 

(k 2 -n 2 ) 2 



(_1)*-H»y(l)( ff ( i ))_y(l)(0) 



(A.7) 



Appendix B. WKB approximation 

In this appendix we describe the application of the WKB approximation to 
the problem of a string of variable density and derive the first few corrections 
to the energies of the modes of the string. Our discussion follows closely ref. Q , 
where however only the leading contribution was derived explicitly. 

Our starting point is eq. (JIJ: we assume 

V n (x) pa e L ^^=o^S n (x) ^ ( B1 ) 
with 5 — > 0. Assuming E n to be of order 1/5 2 , eq. (QJ may be written as 

- 5 2 ^P- = E n p(x)* n (x) . (B.2) 
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Substituting eq. (|B.1|) inside eq. (|B.2|) we obtain the equations for the dif- 
ferent orders in 5. To leading order, for example, the equation obtained is 



whose solution is 



S' {xf + E n p(x) =0 



S (x) = ±i\/E^ J \/ p(y)dy 



To first order we have 



whose solution is 



2S' (x)S[(x)-So(x) =0 



Si(x) = --logp(x) + c 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



where c is a constant of integration, whose value should be fixed by the normal- 
ization of the solution. 

To second order we have 



S[(x) 2 - 2S' {x)S' 2 {x) - S'{{x) = 



(B.7) 



whose solution is 



S 2 (x) = -J 



* S[{yf + S'{{y) 



f x 5p'(y) 2 - My)p"(y) 



-L 

i 



dy 



(B. 



Substituting these results in eq. (|B.lj) we have 

or equivalcntly 

* n (x) cx 



p(x) L > 4 



1 



p(x) 1 /* 

1 



sin 



32 JE. 



E n J ^ Vriyjdy 
* 5p'(yf-4p(y)p"(y) 



n J -L 



P(y) 5/2 



dy 



(B.10) 



where the sign ambiguity is irrelevant and can therefore be removed. 

We need to enforce the boundary condition at a; = L, thus obtaining the 
" quantization" condition 



1 



/ Vp(y)dy- 

<-L 6L\J hi n J-L 



L 5 P '(y) 2 - Mv)p"(y) 
P {yf /2 



dy = nir . (B.ll) 
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Although this equation may be solved exactly it is convenient to assume 

oo 

v/^-^TMmr) 1 -" , (B.12) 



n=0 

and find the constant coefficients a n : 



a = —f J=^- = -jjr , (B.13) 

ai = (B.14) 



Therefore: 



En ~^W + ^L)J_ L 16K^ dy + (BJ6) 

This expression agrees with the one found with WKBPT. 
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